AD-A173  820 FINITE  ELEMENT  ANALVSIS  OF  PROBLEMS  IN 'I 
ELASTIC  I  TV/PLASTIC  I  TV  WITH  SINGULARITIES^)  8RUNEL  UNIV 
UX8RIDGE  i ENGLAND  >  INST  OF  COMPUTATIONAL  MATHEMATICS 
UNCLASSIFIED  J  R  WHITEMAN  SEP  86  DAJA45-8S-C-0855  F/G  2«/il  NL 


\ 


!CROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURfAU  Of  STANDARDS  1%3-A 


AD- A 173  820 


FINITE  ELEMENT  ANALYSIS  OF  PROBLEMS  IN 
ELASTICITY /PLASTICITY  WITH  SINGULARITIES 

by 

J.R.  Whiteman 


September  1986 


DTIC 


United  States  Army 

EUROPEAN  RESEARCH  OFFICE  OF  THE  U.S.  ARMY 
London  England 

CONTRACT  NUMBER  DAJA45-83-C-'t)55  R(D  4238  MA 


Institute  of  Computational  Mathematics,  Brunei  University 


Approved  for  Public  Release;  distribution  unlimited 


FINITE  ELEMENT  ANALYSIS  OF  PROBLEMS  IN  ELAST1C1TY/PLAST2 CITY 


WITH  SINGULARITIES 


1 •  AIMS  AS  STATED  IN  THE  ORIGINAL  PROJECT  PROPOSAL 

The  aims  of  this  project,  as  set  out  in  the  original  proposal,  may  be 
summarised  as  follows:  an  investigation  would  be  undertaken  into  the  use  of 
finite  element  methods  for  three-dimensional  linear  problems  from  potential 
theory  and  elasticity  containing  singularities  and  for  nonlinear  problems, 
particularly  those  involving  the  elasto-plastic  deformation  of  three- 
dimensional  solids.  Special  emphasis  would  be  given  to  the  forms  of 
singularities  in  three-dimensional  linear  problems  and  the  use  which  could  be 
made  of  these  in  finite  element  methods,  to  the  use  of  complementary  energy 
formulations  and  variational  inequalities  in  elasto-plastic  problems  and  to  the 
treatment  of  elasto-pl ast ic  problems  involving  cracks.  Throughout  consider¬ 
ation  would  hr  ,-iven  to  the  theoretical  basis  of  t  lie  finite  element  met !  »-ds 
dev*.  1  oped  for  these  pioble.'.s,  U  tie  development  of  effective  algor  i  t  Lmt  for 
their  imp] omenta t i on  and  application  to  solve  practical  problems. 

Tiicsc  tasks  have  boon  adhered  to  and  the  results  of  the  research  are  outlined 
below.  Both  linear  and  nonlinear  problems  have  been  studied  and  these  arc 
treated  separately  in  Section  2  (linear)  and  Section  3  (nonlinear).  A  total  of 
cl  even  puol  iiat  ions  !  nve  resulted  from  this  re  search .  These:  are  listed  in 
Section  7  and  a:iu  should  be  read  in  conjunction,  with  this  Report.  ^ _ 

2 .  LINEAR  PROBLEMS 

As  was  outlined  in  the  original  proposal,  whilst  many  theoretical  results  exist 
concerning  the  form  of  singularities  in  two-dimensional  problems  of  potential 
theory  and  linear  elastic  fracture,  much  less  is  known  concerning  the  form  of 
singularities  for  three-dimensional  problems  of  this  type.  For  example  in 
opening  mode  linear  elastic  fracture  there  is  still  considerable  discussion 
regarding  the  form  of  the  singularity  at  the  point  where  a  crack  front  meets  a 
stress  free  surface.  Vnen  finite  element  techniques  are  used  to  treat  problems 
involving  singularities  their  numerical  efficiency  and  accuracy  depends  upon 
the  manner  in  which  the  numerical  solution  is  able  to  model  the  singular 
solution.  It  has  long  been  known  that  standard  finite  element  techniques  do 
not  model  singular  solutions  effectively,  so  that  many  special  adaptations  have 
been  made  to  improve  this  situation.  Whatever  special  technique  is  used,  it 
relies  on  some  knowledge  of  the  form  of  the  singularity.  Thus  as  in  the  above 
three-dimens iona 1  contexts,  when  this  form  is  not  known  it  is  not  possible  to 
produce  an  effective  finite  element  adaptation  with  error  analysis  for  the 
singular  problem. 

One  major  feature  of  the  work  of  this  project  has  been  the  theoretical  analysis 
of  three-dimensional  Poisson  problems  containing  re-entrant  vertices.  In 
addition  to  this  the  use  of  singular  forms  in  finite  element  techniques  for 
such  problems  has  been  studied.  On  the  theoretical  side  Beagles  and  Whiteman 
have  produced  two  papers  [1]  and  [3]  which  treat  t hree-d i mens i ona 1  Poisson 
problems  and  derive  both  theoretically  and  numerically  the  forms  of  vertex 
singularities.  For  certain  problems  such  as  those  of  Figure  1  the  singular 
forms  cannot  be  derived  exactly  and  numerical  approximations  have  to  be 
obtained.  In  til  and  1 3 1  numerical  schemes  are  given  for  obtaining  such 
approximations  by  solving  Lap  1 acc-Bel t rami  problems  and  eigenvalues,  which  give 
t  lie  exponent  of  the  radial  behaviour,  have  been  calculated  together  with  upper 
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and  lower  bounds.  When  these  calculated  values  have  to  be  used,  the  singular 
functions  are  "non-exact"  and  so  study  of  the  use  of  non-exact  singular 
functions  of  this  type  in  finite  element  techniques  has  also  been  made, 
see  [2],  As  a  result  of  all  this  work  we  now  have  a  much  better  understanding 
of  the  forms  of  vertex  singularities  in  three-dimensional  Poisson  problems  and 
their  use  with  finite  element  techniques.  A  wide  range  of  vertices  has  been 
treated  in  [3] . 

In  the  context  of  linear  elasticity  Whiteman  and  Thompson  have  studied  the 
stress  forms  for  three-dimensional  problems  containing  cracks  and  notches  and 
have  obtained  some  preliminary  results  concerning  such  forms.  This  work  has 
again  involved  the  solution  of  eigenvalue  problems  and  a  manuscript  til]  has 
been  written.  As  this  subject  is  one  of  particular  importance  and  is  somewhat 
an  area  of  dispute,  copies  of  this  manuscript  are  being  sent  to  various 
colleagues  who  work  in  this  field,  prior  to  the  manuscript  being  submitted  for 
publication.  In  summary  there  are  two  "camps"  in  this  subject,  one  of  which 
claims  a  ccitui.i  i..i  ..  : .  1  Liu..  singularity  and  one  of  which  claims  another.  An\ 
»«  w  paper  on  t  Hr  suhj.-it  ic  thus  bound  to  be  slightly  cent  rover:-  i  -3 1  .  Or.cc  this 
r.anusc  r  i  t  is  finished  it  i  my  intention  to  int  orpot  ale  our  singular  fields 
into  our  finite  clement  codes  and  to  undertake  analyses  similar  to  those  which 
have  been  performed  for  the;  three-dimensional  Poisson  problems. 

In  the  context  of  two-  and  three-dimensional  linear  elastic  fracture  Thc.'pson 
and  1  tern:.  *.  -ins 1  acted  lie  family  oi  quadratic  ele:i;its  w.ich  arc  modi  i  i  cd  to 
produce  singular  approximating  forms  by  the  repositioning  of  certain  nodes; 
usually  on  the  element  edges  which  intersect  at  the  point  of  singularity.  This 
is  the  family  of  "quat ter  point  elements"  and  includes  transition  elements. 

In  [4]  a  detailed  analysis  is  given  of  the  strain  forms  produced  with  these 
elements  at  and  near  crack  tips.  This  indicates  the  strain  forms  for  various 
well  known  elements,  particularly  radial  strains  along  lines  emanating  from  the 
point  of  singularity.  The  paper  [4]  sets  out  to  show  when  the  correct  singular 
form  is  produced  and  when  it  is  not.  It  has  since  been  used  by  Michavila  and 
his  co-workers  in  Madrid  [14]  as  a  prototype  for  the  analysis  of  other  singular 
el ements  . 

Superc  onvcrgence 

An  aspect  of  the  finite  element  treatment  of  linear  problems,  which  was  not 
envisaged  when  the  proposal  for  this  contract  W3S  written  and  which  I13S  emerged 
during  the  contract  period,  is  that  of  super c onvcrgence  in  finite  element 
techn  iques  . 

It  has  been  known  for  some  years  that  for  finite  element  solutions  to  Poisson 
problems  and  problems  of  linear  elasticity  there  are  specific  points  (stress 
points)  at  which  accuracy  is  better  than  that  found  globally.  In  t he  last 
three  years  several  papers  have  appeared  which  give  theoretical  analysis  and 
proofs  that  gradients  can  be  recovered  from  finite  element  solutions  to  Poisson 
problems  in  such  a  way  that  at  specific  points  the  rate  of  convergence  is 
higher  than  that  found  overall.  This  is  the  phenomenon  of  superc onvcrgence . 
Unfortunately  the  main  shortcoming  of  these  superconvergence  estimates  is  that 
they  require  high  global  regularity  of  t be  solution  of  the  original  problem; 
exactly  the  property  that  the  solution  of  a  boundary  value  problem  with  a 
singularity  does  not  possess. 
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In  association  with  Professor  M.F.  Wheeler  of  Kice  University,  I  have  devoted 
considerable  time  to  deriving  error  estimates  for  the  superconvergent  recovery 
of  gradients  for  finite  element  approximations  to  Poisson  problems  in  two 
dimensions.  Our  major  contribution,  which  appears  in  [5]  concerns  the  super¬ 
convergence  on  sub-domains  of  the  problems.  By  deriving  results  on 
sub-domains,  where  the  solution  has  high  local  regularity,  we  are  able  to  treat 
problems  involving  singularities.  We  consider  the  results  of  [5]  to  be  a  very 
significant  advance. 

The  superconvergence  results  requiring  high  global  regularity  for  Poisson 
problems,  which  have  been  published  previously  in  the  literature,  have  now  been 
adapted  by  Whiteman  and  Goodsell  to  the  case  of  linear  elasticity  and  it  is 
intended  to  adapt  these  in  turn  to  the  case  of  linear  elastic  problems  in¬ 
volving  singularities  such  as  occur  in  linear  elastic  fractuie.  Particularly 
significant  will  be  the  ability  to  obtain  superconverger.ee  results  for  cal¬ 
culated  values  of  J- integral  s  .  T  n  view  of  the  widespread  u-e  of  J- i  u  t  «-g  i  u  1 1  as 
fracture  criteria,  such  an  advance  would  be  very  important . 

_  3.  NON'Ll  .\EAK  PROBLEMS 

The  second  major  feature  of  the  project  has  been  the  treatment  of  problems  of 
e 1  a h t o-p 1  a s t i c i t y .  In  tin-  proposal  a  formulation  of  an  el  a st it -plastic  problem 
in  t« rms  of  a  variational  inequality  by  Oden  and  Whiteman  was  set  out  [12]. 

This  was  formulated  in  terms  of  complementary  energy  and  was  based  on  the 
deformation  theory  of  plasticity.  One  stated  possibility  was  to  adapt  this  to 
treat  problems  involving  elastic-plastic  fracture.  Inspite  of  much  effort  and 
collaboration  between  Oden  and  Whiteman,  the  complementary  energy  formulation 
in  terms  of  stresses  proved  unsuitable  for  the  fracture  context.  From  the 
engineering  point  of  view  this  could  perhaps  have  been  expected  in  view  of  the 
shortcomings  of  deformation  theory  plasticity.  The  difficulty  stems  from  the 
fact  that  stresses  are  the  primary  variables,  and  for  the  linear  elastic 
fracture  model  these  are  infinite  at  a  crack  tip. 

In  order  to  treat  e 1  as t o-pl as t ic  fracture  problems,  we  therefore  considered 
incremental  plasticity  and  adapted  the  MODEL  finite  element  code  to  treat 
problems  of  nonlinear  fracture  mechanics,  see  [8].  This  code  based  on  the 
incremental  plasticity  model  works  with  displacements  as  primary  variables  and 
in  the  linear  elastic  context  uses  a  potential  energy  formulation.  Once  the 
code  produced  reasonable  results  for  displacements  in  a  number  of  problems,  it 
was  extended  to  calculate  integrals  of  the  J-type  appropriate  to  the  context  of 
elastic-plastic  deformation  for  problems  with  cracks.  Several  papers,  [b  -  1 0 ] 
were  written  on  this  aspect  and  calculated  values  of  this  path  integral  were 
derived.  It  was  found  from  the  numerical  results  that  for  the  nonlinear  case, 
where  significant  elastic-plastic  deformation  takes  place,  the  specific 
integral  appears  to  be  no  longer  path  independent.  During  a  visit  to  MTL, 
Watertown,  in  May  1983  I  had  lengthy  discussions  with  Dr.  D.M.  Tracey  who 
confirmed  his  belief  that  this  effect  is  correct.  Work  is  now  underway  to 
incorporate  other  integrals  of  the  J-lype  into  the  code.  These  last  integrals 
were  proposed  in  [ 1 5 ]  and  are  path  independent  in  this  context. 

Thus  the  current  situation  is  that  the  MODEL  code  has  been  extended  to  treat 
nonlinear  fracture  for  problems  involving  elastic-plastic  deformation  and  has 
been  tested  on  two-dimensional  problems.  It  lias  also  been  applied  to  a  limited 
number  of  t hree-d imen s i ona 1  problems.  In  this  latter  case,  since  integrals  of 
the  J-typc  are  two-dimensional  phenomena,  a  family  of  J  values  has  been 


5 


calculated  for  planes  in  the  solids  orthogonal  to  the  crack  front.  Future  work 
will  involve  the  seeking  of  a  J-type  quantity  appropriate  to  nonlinear 
three-dimensional  problems  of  this  type.  All  this  has  been  done  by  numerical 
experiment.  At  the  same  time  I  have  been  seeking  to  formulate  a  rigorous 
mathematical  model  of  the  incremental  plasticity  problem  with  hardening,  which 
can  be  used  to  provide  a  theoretical  basis  with  error  estimates  to  the  finite 
element  method  in  this  context.  Using  the  analysis  of  Korneev  and  Langer  [ 1 3  3 
1  now  have  one  such  model.  Discussion  with  Professor  J.T.  Oden  has  taken  place 
over  this.  The  next  step  is  to  perform  the  finite  element  analysis  for 
elastic-plastic  problems  without  singularities  and  then  to  extend  this  to  the 
nonlinear  fracture  case. 


4.  CONCLUSIONS 

The  work  of  this  project  may  be  summarised  as  follows: 

(a)  A  considerable  number  of  results  have  been  derived  concerning  the  forms  of 
singularities  in  three-dimensional  Poisson  problems. 

(b)  Very  significant  results  have  been  derived  for  the  superconvergent 
recovery  of  gradients  of  finite  element  solutions  to  Poisson  problems.  It  is 
anticipated  that  these  will  he  important  for  use  in  fracture  mechanics  in  the 
context  of  the  calculation  of  path  integrals;  immediately  for  linear  p rob It. as 
but  ultimately  in  nonlinear  fracture  mechanics. 

(c)  Software  now  exists  for  the  treatment  of  nonlinear  fracture  problems  and 
it  can  be  used  to  calculate  J-type  integrals  in  this  context.  A  theoretical 
model  and  finite  element  error  estimates  are  being  derived.  The  next  stage 
will  be  to  derive  and  compute  an  integral  of  the  J-type  for  use  on  cracks  in 
ihree-dimcnsiona 1  problems  involving  elastic-plastic  deformation. 

(d)  The  collaboration  with  Professor  J.T.  Oden  of  the  University  of  Texas 
proved  very  useful  from  the  point  of  view  of  discussions  concerning  the 
treatment  of  e 1  a s t i c/pl a s t i c  problems  in  the  context  of  mathematical  formul¬ 
ations  for  treating  plastic  behaviour  of  this  type.  Extensions  to  the  non¬ 
linear  fracture  case  are  now  being  considered. 

5.  FUTURE  WORK 


As  is  usual  at  the  end  of  any  project,  many  questions  remain  unanswered  and 
many  new  ones  have  been  posed.  Clearly  exciting  opportunities  exist  for  the 
exploitation  of  the  local  superconvergence  results.  The  calculation  of 
J-quantities  which  are  path  independent  in  the  nonlinear  context  and  the 
derivation  of  an  integral  for  three-dimensional  nonlinear  fracture  are  also 
extremely  important.  During  the  course  of  the  project  my  attention  has  turned 
to  more  complicated  materials  such  as  polymers,  where  phenomena  which  are 
tine-dependent  come  in  and  also  where  non-isotropic  effects  are  present.  1  am 
tiius  currently  studying  the  fracture  of  orlhotropic  materials  and  also  the 
fracture  of  materials  where  viscoelastic  and  viscoplastic  deformation  occurs. 

In  this  last  context  1  believe  it  is  significant  that  1  have  been  invited  to 
spend  the  Fall  Semester  of  1986  as  a  Visiting  Professor  in  the  Department  of 
Mathematics  at  Texas  A  &  M  University,  where  I  shall  collaborate  with 
Dr.  J.  Walton,  and  also  with  Professor  R.  Shapery,  who  is  the  Director  of  the 
Mechanics  and  Materials  Center,  and  who  is  in  addition  a  world  authority  on  the 
theory  of  viscoelastic  fracture. 
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6.  LECTURES  AND  VISITS 


During  the  period  of  this  contract  I  was  invited  to  present  talks  at  the  3rd 
and  4th  US  Army  Conferences  on  Applied  Mathematics  and  Computing  and  1  am  most 
grateful  to  the  US  Army  European  Research  Office  in  London  for  their  support 
which  enabled  me  to  attend  these  conferences.  In  addition,  during  the  period 
September  1983  to  September  1986,  I  have  visited  the  Benet  Laboratory, 
Watervliet,  and  the  MTL,  Watertown,  on  several  occasions  and  on  each  visit  have 
had  very  useful  discussions.  A  list  of  the  main  lectures  which  I  have 
jj  presented  in  various  parts  of  the  world  during  the  contract  period  is  given  in 

Appendix  I. 
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APPENDIX  ] 


Main  Lectures,  September  1983  -  September  1986 


GAMM  Workshop  on  Numerical  Analysis,  Kiel 

MAFELAP  1934,  Conference  on  The  Mathematics  of  Finite  Elements 
and  Appl icat ions ,  Brunei  University,  Uxbridge 
Summer  School  of  Numerical  Analysis,  Lancaster 
Tochnische  Hochschule,  Hannover 
Techmsche  Hoc h schu  1  e  ,  Darmstadt 


Conference  on  Mathematical  Ar.alvei  s  and  Applications,  Kuwait 
Ki  net  J  ;,-hor  3  t  orv  ,  U'afprv'  iaf 

D*  j  *  a  r  t  me  n  t  ■  >  i  C’-’  ■  t  c  r  Sc  auk  r  ,  ho  n  s  s  o  ]  a  e  r  Polytechnic 

1:  1  ;  ••  A:-.y  .'-ppi  it  ..  Me  t  S.  Co- put  i-.c  Ler.f  •:  r<  ;.ce  ,  At  1  n*.t 

Department  of  Mat  j.e.rat  its  ,  University  of  T->...s  at  A.ist  i  n 

Di-par  f  one  of  Mathematical  Sciences,  Kit  e  University 

Department  of  Mathematics,  University  of  Wales,  Aberystwyth 

Conference  on  Numerical  Analysis,  Dundee 

Conference  on  Numerical  Analysis,  Madrid 

Workshop  on  Discretisation  in  Continuum  Mechanics,  Bad  Honnef 


Material  Technology  Laboratory,  Watertown 
Department  of  Mathematics,  Texas  A  4  M  University 
Department  of  Mathematical  Sciences,  Rice  University 
Department  of  Mathetnat  ics ,  Herriot  Watt  University 
Edinburgh  Mathematical  Society 

Department  of  Mathematics,  University  of  Kuwait 
Department  of  Mathematics,  University  College  of  Bahrain 
School  of  Mather, at  ics ,  University  of  Bach 

Fourth  Army  Applied  Mathematics  4  Computing  Conference,  Atlanta 

Department  of  Mat henat i c s ,  University  of  Strathclyde 

School  of  Mining  Engineering,  University  Polytechnic  of  Madrid 


January  1934 

May  1984 

July  1934 

September  1984 
St  p  1 1  mb*  r  3  9S4 


February  19  85 
May  1985 

•‘•ay  1  Pop 


X  ay  1  9  o  p 

May  1985 

June  1985 

September  1985 
September  1 985 


J a nua  ry 

1986 

January 

1986 

Jam. a  ry 

1986 

March 

1986 

March 

1986 

April 

1986 

Apr  i  1 

1986 

May 

1986 

May 

1986 

June 

1986 

July 

1936 
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